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1 Introduction 



Over the past decade information theory has been generalized to include quantum mechanical 
systems, for example, a two-level quantum system has come to be known as a qubit in this context. 
The additional freedom introduced with the quantum mechanical superposition principle has opened 
up a variety of capabilities that go well beyond those of conventional information techniques. There 
are two distinct directions in which progress is currently being made: quantum computation and 
error correction or prevention on the one hand[l], and nonlocality and distillation, on the other 
hand [2]. In each of those progresses, quantum entanglement that provides a good measure of 
quantum correlations plays an important role. 

There are a number of good measures of the amount of entanglement for two quantum systems 
in a pure state, a good measure of entanglement for mixed states is also found though it is hard to 
compute for a general state [3]. 

In present paper, considering the entanglement from the other aspect, we prefer to discuss the 
change of entanglement due to the state changes rather than to compute straightforwardly the 
entanglement of an arbitrary state. 

In the framework of quantum information theory, the state change allowed by quantum mechan- 
ics may by treated in terms of quantum operations [3], a simple example is the unitary evolution 
experienced by a close quantum system. The final state of the system is related to the initial state 
by a unitary transformation U, 

p^e(p) = U P U + . 

Unitary evolution is widely in use of quantum gates and circuits [4] as a quantum operation. In 
addition to the unitary evolution, environment coupling to a quantum system or a measurement 
performed on the quantum system changes the state too[6,7]. The connection of quantum operations 
to quantum measurements is easy to explain. Standard text book treatments describe quantum 
measurement in terms of a complete set of orthogonal projection operators for the system being 
measured. This formalism, however, does not describe many of the measurements that can be 
performed on a quantum system. The most general type of measurement that can be performed 
on a quantum system is known as generalized measurement [6] . Generalized measurement can be 
understood within the framework of unitary evolution, because most generalized measurements can 
be realized through many dynamical processes [8], and the state change due to environment may be 
also treated in terms of the unitary evolution, since an arbitrary open system may be enlarged by 
including the environment to be a close system. In this sense, unitary evolution is one of the most 
general types of state change possible in quantum mechanics. 

The rest of present paper is organized as follows: In Sec. 2, we present a general method to fac- 
torize the unitary evolution operator U (t) for a close system. The results may be generalized in the 
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treatment of many autonomous dynamical systems. Sec. 3 contains our results on the entanglement 
change occurring in a dynamical process. Finally, in Sec. 4, we present two typical examples and 
some conclusions. 



2 Factorizing the unitary evolution operator U(t) 

As noted above, the unitary evolution operator is one of the most general types of state change 
possible in quantum mechanics, the point of this section is to factorize the evolution operator into 
a set of independent one. For this end, we discuss the following cases. 
case A [9] — The Hamiltonian can be written as a finite sum 

m 

H{t) = Y,<t)Hi, (2.1) 

i 

where Oj(t) are a set of linearly independent complex valued functions of time, and Hi are constant 
operators. In addition, the set of operators Hi(i = 1, ...,m) may be enlarged by repeated commu- 
tation to a Lie algebra L with finite dimension n(n > m). With this presupposition, the unitary 
evolution operator can be uncoupled into a set of independent operators 

U(t) = U 1 (t)U 2 {t)...U n {t), (2.2) 

where each component Ui(t) is an operator satisfying 

j t Ui(t) = gMHiUiit), Ui(0) = 1. (2.3) 

With the scalar function gi(t) being the solution to a set of nonlinear differential equations 

d 11 

j t 9i(t) = J2 r lika k (t),g i (0) = 0, (2.4) 

i=l 

where 77^ are nonlinear function of g's. Thus we have factorized the unitary evolution operator 
into the form: 

n 

U(t) = '[[e!> i ® Hi . (2.5) 

i=i 

Especially, for a general case of a dynamically closed quantum system which consists of two inter- 
acting subsystems A and B, the total Hamiltonian may be written as a sum of three terms 

H = H A + H B + H int , (2.6) 

the first two terms represent the free Hamiltonian of subsystem A and B, respectively, and the last 
term describes the interaction between the A and B. Following the procedure stated above, we 
arrive at 

U(t)= J] e 9x{t) Hi. (2.7) 

i=AB,int,...,M 
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Here, Hi(i = A, B, int, ...M) are elements of Lie algebra with finite dimension enlarged by Ha, HB,Hi nt 



case B — In the case of the dimension of the Lie algebra enlarged by {H{\ is infinite, we can fac- 
torized the unitary evolution operator using the general Baker-campbell-Hausdorff[10] formula. To 
start with, we give the evolution operator of the system under consideration 

U(t) = e~ itH = e -it(HA+H B +H int ) ^ ^ 

where Ha, Hb, Hi nt are the same as in eq.(2.6), the eq.(2.8) can approximately be written as 

U{t) = e -^l H A+H B ,H mt ]t 2 . e -iH int t . e -i(H A +H B )t _|_ 0(t 3 ). (2.9) 

This splitting formula is hold in the case that t has to be safely smaller than a typical energy of the 
system. Thus, even in the simplest case, a better method is needed. Let n be a positive integer. 
The exponential function satisfies the scaling identity 

exp(-iHt) = [exp(-iHt/2 n )] 2n . (2.10) 

When n is sufficiently large, the argument ^ is in some sense small. Eqs.(2.9) and (2.10) together 
give 

JJ(t) = e -k[ H A+H B ,H lnt \T 2 . e -iH int r . e -i(H A +H B )T e -\[H A +H B ,H in t\r 2 . 

e -iH lnt r . e -i(H A +H B )r ^^ e -±{H A +H B ,H int ]T 2 . e ~iH int T . e ~i(H A +H B )r _j_ 0(t 3 ), (2.11) 

where r = Still higher-order formulae are known. We would like to point out that, in quantum 
computation [11], the n = 1 is widely taken in use and it is large enough to avoid the decoherence 
during quantum computing. 

3 Quantification of entanglement 

In the previous section we have factorized the time evolution operator U(t). The question remains 
open about how does the entanglement change in a dynamical process. Of course, this question is 
not entirely well defined unless we state what physical circumstances characterized the amount of 
entanglement. This suggest that there is no unique measure of entanglement. Before we define the 
measure of entanglement we expand the density operator for a close system. Suppose that the two 
interacting subsystems are initially separable[12], i.e., their initial density operator (state) can be 
written in a form 

p(0) = p A (0) ® p B (0), (3.1) 

to use the entanglement for quantum information processing, however, we need a inseparable state, 
more precisely, a state in pure entanglement form. The procedure of converting a separable state to 
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inseparable one can be performed, as stated in Sec.l, through a unitary evolution operator J7(i)(in 
addition, a partial trace is also needed sometimes.) 

p(t) = U(t)p(0)U + (t). (3.2) 

If the interaction between the two subsystems is small, it is natural to attempt some sort of Taylor 
series expansion of the exponential in eqs. (3.2) (2.8) and (2.11), which give 
case A 

p(X,t) = p° A (t) ® p° B (t) + \ ]T MHMt) <g> (P B {t) + (P A {t) ® p%{t)^H$ 

i^A,B 

+ y E (^^H t HA(t)® PB (t) + p A (t)®p B (t)^^H iHj 

+ ^H t p A {t)®p%{t) d -^-H 3 ) + 0(\ 2 ) (3.3) 



and case B 

p{\t) = p A (t)®p° B (t)-^) 2 ^{[H A + H B ,H int (t i )Up A (t)®p° B (t)}+ 

i=0 

2 n — l 

- * A (^) Y,i H ™t(U),PA(t)®P B (t)]-+0(\ 2 ) + 0((^f), (3.4) 

i=0 

where A denotes the coupling constant, and pf(t) represents the state of subsystem i at time t with 
A = 0. The results presented above suggest that we may take the form 

5D(p) = \\p{t) - p° A (t) ® p° B (t)\\ 2 = Tr(p(t) - p° A (t) p%{t)f (3.5) 

as a measure of entanglement change in the time evolution process. Noticing the initial state 
is separable, the measure of entanglement change (3.5) is a measure of entanglement in reality. 
Although the definition of the measure for entanglement is not unique, it has to satisfies the three 
condition stated below [3]: 

(1) .D(p) = if and only if p is separable. 

(2) Local unitary operators leave D(p) invariant, i.e. 

D(p) = D(U A ® U B pU+ ® U+). 

(3) The expected entanglement cannot increase under Local general measurements+ Classical com- 
munication+Postselection (LGM+CC+PS) given 

by EiV+Vi, = l,i.e., 

Y / Tr(p l )D(^/Tr Pt )<D(p), 

wherepi = VipV^ . For the measure of entanglement change proposed above, (l)follows from the 
fact that D{p) is a true metric, and (2) is obvious. Property (3) is satisfied too [3]. We believe that 
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there are numerous other nontrivial choices for measure of entanglement, one of the choices could 
not be said to be more important than any other,the present choice has the advantage that it is 
easy to compute for any dynamical process. 

Our discussion so far has centered on the entanglement change in a dynamical process. To 
complete it we still need to show that this definition can be generalized for any process that 
quantum mechanics allowed. For a general process, the quantum operator that change a state of 
the system should be factorized by the subsystem's operators. For instance, a control not operation 
in quantum computation given by (in fact, control not is a unitary evolution operator) 



where I2 is the unit operator for the second qubit, erf stands for the x pauli matrix of the second 
qubit. I l)i and |0)i represent two state of the first qubit. This control not operator is factorisable, 
i.e. O can be written in the form 



where 0\ and 0\ denote operators for the first and second qubit, respectively. Hence, according to 
the definition eq.(3.5), the control not operator in the form (3.6) does not change the entanglement 
of the system. 

We would like to point out that the discussions presented here are for the unitary evolutions, 
for non-unitary evolution such as a trace over some of the degree of freedom, we should find a 
auxiliary unitary process instead of the non-unitary one. 

4 Example and Conclusion 

In order to understand how our program for calculating the amount of entanglement change works, 
we present in this section two examples, one of them consists of two interacting qubits (two identical 
two-level system) in a laser beam [5a], and the another two independent qubits coupling simultane- 
ously to a bath. 
example 1: 

The Hamiltonian describing the system in this example has the following form(set U = 1): 



O=|0)i(0|®l 2 + |l>i<l|®af, 



(3.6) 



= E°i®°2 



H = Ha + Hb + Hi n t + Hf, 
Hi = -u<ri,(i = A, B) 

Hint = g(<?\<?B + u ~a°b) + A ( S 



aja + a ia ), 




(4.1) 



6 



where of ,0^ ,af describe the pauli operator of the i qubit, g denotes the coupling constant, and 
Hf stands for the free Hamiltonian of the laser beam. Suppose the state is initially in the form 

p(0) = P A(0)®p B (0)®Pf(0), 
Pa(0) ® Pb(0) = \e A ,e B )(e A ,e B \, 

p/(o) = Epwwh, ( 4 - 2 ) 

n 

where |ej) denotes the excited state of the qubit i and \n) stands for a Fock state of the laser beam 
field. In the Schrodinger picture, the density operator that obeys von Neumann equation is given 
by 

p AB (t)=Tr f p(t) = 

( S„ p2(n)/ 9S (n ' t) D„ p(n+ 1 )P(")/ss( n + 1>*)/EG(™,*) XL P(U+ 2 )P(™)/ee(n+ 2)/ 3s (n) \ 

S„ p(n + 1)p(n)/ " (n + 1 '* )/EG(n '* ) S„ p2 ( n)/ EG( n '*) ^p(n)p(n + l)/ EG (n+ l)/ee(») , 

V X)„ p(n + 2)p(n)/ " (n + 2)/ 99 (n) X)„ p(n)p(n + 1)/EG(n + 1)/ " (n) X)„ p2(n)/ " (n,t) / 



(4.3) 



where we take |5a,5s), G), and \eA,e B ) as a set of basis, and 

\9a,9b) = \9a) ® \9b), \eA,e B ) = \e A ) <8> |e B ), G) = ^(1^) ® |e B ) + (e^) 

P^)®Pb(*) = Pa(0)®p§,(0), 
/ ffS (n, i) = J sin 20 sin fle - ^' + i sin 2<j) cos 2 ^e~ iB -* - i sin 2^ e - iEot , 

f ee (n, t) = sin 2 sin 2 ^e~ iS +* + cos 2 4>e~ iE ^ + sin 2 cos 2 \z~ iE ~\ 

fEG(n,t) = sin sin sin 

_, „ cos ± 1 , _, . ,_, . _, . 
£± = n + w/(n + 1), #o = (n + IV/, 



and0 = § ,tan^> = y'^if ,17 2 = (16n + 24)c/ 2 . Eqs.(3.5) and (4.3) together give 

5D(p)= E (pSb) 2 " + 1, (4.4) 

ij=l,2,3 

where p %3 AB denotes the element of matrix pa B given by eq.(4.3), which represents the entanglement 
change or entanglement of subsystems A and B at time t. 
example 2 

The Hamiltonian describing dissipation of the two qubits has the following formfld] (setting U = 1) 
H = u ((j z a + a z h ) + E / M9u,lMrti + Owl)] 

+ y da; U J=0)6 (wa+ja^), (4.5) 
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where &i describe the pauli's matrix of the i qubit, a w i stands for the bath mode to coupling to 
the I qubit, and yJi= a ,baZl a ul = a tl a ^l for «™ = a ub , whereas Uz =a , 6 a^a^ = a+ a a ua + a^a^ for 
a-Laa 7^ ciuib- The coupling coefficients are denoted by g w i, and the qubit operator A\ in general is 
expressed as a linear superposition of three pauli's operators, i.e. A\ = A^Vf + A^of + 

A (3) 

or. 

The ratio A^ : A^ 2 ) : A^ 3 ) is determined by the type of the dissipation. For instance, A^ = A^ 2 ) = 
for phase damping and A® = for amplitude damping[7]. Phase damping induces pure dephasing, 
whereas amplitude damping induces loss and dephasing simultaneously. Many source of decoherence 
in quantum computers are described by amplitude damping[13]. 

Without any loss of generality, we discuss in detail the case with A^ = A^ 2 ^ = i.e. phase 
damping. Some words of caution are now in order. As mentioned above, the bath can also cause 
another unwanted effect in computation process, i.e. amplitude dissipation. It is easy, however, 
to make system have small loss rate of amplitude dissipation [11], so a considerable number of 
operations are allowed to perform. 

In the case of phase damping, the unitary evolution operator may be factorized in the following 
form 

U(t) = U (t)Uj(t), (4.6) 

where Uo(t) = e~ tHot with Hq = loq{ct^ + a'() + J du> yJi =a ,b {^( 1 t>i a ^i) ls ^ ne ^ ree evolution operator, 
while Ui(t) denotes the evolution operator in the interaction picture. A readily calculation shows 
that 

U I (t) = Uf(t)U b I (t), 

Uj(t) = >< a\ + v\\gi >< gt\, (4.7) 

where |ej > and \gi > are the eigenstates of af with eigenvalues +1 and —1, respectively, and u\ 
and v\ satisfy (% = a,b): 

i% dl ui l = Y, duJ 9coi{a^e iujt +a L0i e- iujt )v i I , 

lh §-/i = -Y. d ^K t e^ + a M e-^)v\. (4.8) 

The Wei-Norman's algebraic method [9] that provides a way to factorize the evolution operator 
gives 

u * = Y[ e ft(t) e Ai(t)a+ e Bi(t) au , (49) 

and 
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Here, 

Al(t) = -^(e wt -l), 
ijj 

Bt(t) = -(Ai(t)r, 

2 2 

fUt) = -i—t+H(l-e- iu}t ), 

ct(t) = -Ai(t),DKt) = -Bi(t)X(t) = m- 

Now, we turn our attention to compute the reduced density operator of the two-qubit system, first 
of all, we calculate the total density operator, which follows straightforwardly from eq.(4.6) 

p(t) = U o U lP (0)U+U+, (4.11) 

where p(0) denotes the initial density operator (state), which may be written in a separable form 

p(0) = p a (0)®p b (0)®p B (0). 

Here, pi(0)(i = a,b) represents the initial state of qubit i, and p B (0) stands for the initial state of 
the bath. In following, we use the notation, \e a ,gb) to indicate the eigenstates of a z a and erf with 
eigenvalues 1 and —1, Tr B indicate a trace over the bath, and Pi(t) to represent the free two-qubit 
state i.e. 

p°(t)=Tr B U (t)p(0)U+(t). 

With this notation, in a subspace spanned by {|11) = \e a ,eb), |12) = \e a ,gb), 1 21) = \g a ,b2), 1 22) = 
Ida, 9b)}, the state of the two-qubit system at time t takes the form 

/ Pllll P1112 P1121 P1122 \ 

P12U P1212 P1221 P1222 

P2111 P2112 P2121 P2122 

\ P2211 P2212 P2221 P2222 / 

where p ijM = P° ijkl -F ijk i{i, j, k, I = 1, 2), p° ijkl = Tr B (ij\p°(t)\kl), and F ijkl = Tr B T, 2 c ,d,e,f=i( u i)i 
(U^cd) (eflttU^+UiUfj+Ukl). Here, (Uf) t = a (i\Uf\i) a , and \2} a = |e ), |l) a = \g a ). In or- 
der to get more information about the reduced density operator, we make some discussion on the 
quantity F^ k \. It can be easily verified that F^i = 1 for i = k, j = I, while Fijki = F^ u - fovi^k 
and j / I. Moreover, the quantity results from the interaction between the two-qubit system and 
the bath, hence it depends on the states of the bath. Although different bathes result in different 
results Fijki, the physical results discussed here do not rely on the bath. In this sense, we may 
consider a simple case with zero temperature. In this case, F^ki is given that 

F^i = Fijki(t) = e-ri^A+A^pHd^ (4 13) 

where Ay (u, t) = 2 ^^^-^ o.5ut ^ stands for the bath spectrum distribution. The eq.(4.13) 
suggests that F^ki approaches zero with the passage of time except some moments at which 



Pab(t) = Tr B p{t) = 



(4.12) 
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j^[\j{uj,t)/S^ l {uj,t)\p{uj)duj = 0. This attractive results might be used in preventing informa- 
tion loss stored in quantum states. Now we come back to the entanglement change, eqs. (3.5) and 
(4.12) together give 

2 

6D(p ab (t)\\p° a (t)®p° b (t)) = ]T (Pijki-plkiKPkiij-pluj)- (4-14) 

i,j,k,l=l 

In summary, we propose a new method to compute the entanglement change in a dynamical 
process.We see the above treatment in Sec. 2 and Sec. 3 does not refer to specific entangled systems. 
This is a desired property as it makes our measure of entanglement universal. Especially, the results 
yielded by present paper can be easily generalized to more than two subsystems, this is just the case 
of many qubits interacting simultaneously with environment. In addition to the measure stated 
above, the quantum relative entropy defined as 

D(p(t)\\p° A (t) p%{t)) = Tr[p(t)(lnp(t) - lnp° A (t) • p° B (t))} 

and the Bures metric given by 

D(p(t)\\p A (t) p%(t)) = 2 - 2yfF(p,p A ®p%), 

with F(p, p A ®p%) = [Tr[\J p° B <8> P A P\J P°b ® Pa]^] 2 are other good measures of entanglement. With 
this modified definitions, the measures of entanglement can be given in easy way. 
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